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1  .  INTRODUCTION1 

In  performing  parametric  design  studies  on  underwater  bodies  where  the 
flow  properties  of  the  boundarv  layer  are  required,  a  fast-running  procedure 
to  calculate  those  properties  is  a  must.  Another  requirement  of  the  boundary- 
layer  calculation  procedure  is  reasonable  accuracy.  Finite  difference 
procedures  can  be  ruled  out  because  they  require  too  much  computer  time. 

Thus  the  choice  must  be  made  from  the  numerous  fast  running  but  less 
accurate  approximate  methods  which  are  based  on  simplified  forms  of  the 
boundarv-layer  equations  such  as  integral  methods  or  local  similarity 
methods,  c.f.  Chapter  6  of  Rosenhead  [1]. 

The  particular  aspect  of  the  boundarv-laver  problem  addressed  here  is  the 
determination  of  the  suction  distribution  required  to  maintain  laminar, 
attached  flow  in  the  axisvmmetric ,  incompressible  case.  .-As  a  rapid  and 
rasonably  accurate  means  of  solving  this  problem  the  local  similarity  approach 
was  chosen.  Walz  in  1941  [2]  was  the  first  to  use  this  approach  as  a  means  of 
solving  the  integral  form  of  the  momentum  equation.  In  this  method  the 
velocity  profile  is  represented  by  a  self-similar  wedge  flow  profile  where  the 
similarity  variables  are  related  to  the  general  external  flow  variables.  The 
advantage  of  this  method  is  that  it  can  be  used  in  regions  of  very  rapid 
pressure  drop  where  the  methods  based  on  polynomial  profiles,  such  as  that  of 
Pohlhausen,  break  down.  In  1956  Smith  [3]  advanced  a  method  very  similar  to 
that  of  Walz  based  on  piecewise  wedge  flow  solutions  over  finite  intervals  in 
body  arc  length.  Smith's  method  becomes  identical  to  Walz's  method  as  the  arc 
length  shrinks  to  zero.  Smith  applied  his  procedure  to  a  variety  of  two- 
dimensional  and  axisvmmetric  cases  without  suction  and  found  the  accuracy  to 
be  comparable  to  any  of  the  integral  methods.  His  method  has  the  advantage 
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of  being  very  simple  to  apply  and  can  easily  be  extended  to  boundary  layers 
with  suction.  This  report  presents  the  extension  of  Smith’s  method  to  the 
particular  case  where  the  amount  of  suction  is  determined  by  the  critical 
Reynolds  number  (neutral  stability)  curve. 


2.  ANALYSTS 


2. 1  The  Method  of  Smith 

Smith's  method  rests  on  the  assumption  that  a  general  incompressible 
laminar  boundary  layer  solution  can  be  closely  approximated  by  a  piecewise 
matching  of  Falkner-Skan  similarity  solutions  for  wedge  flows.  The  analysis 
starts  with  the  wedge  flow  surface  speed  expression 

3 

u=k(x-x)2‘,  (1) 

e  o 

where  ue  is  the  inviscid  surface  speed,  x0  is  the  arbitrary  origin,  k  is  a 
constant,  8  is  the  Falkner-Skan  pressure  gradient  parameter,  or  Hartree 
parameter  and  x  in  the  present  context  is  body  arc  length. 

Equation  (1)  is  matched  to  the  actual  surface  speeds  at  Points  l  and  2  of 
an  arbitrary  interval  x^  <  x  <  X2<  This  matching  yields  the  two  relations: 

3 

u  =  k(x.  -  x  )2-fl  ,  (2) 

e!  1  0 

B 

u  =  k(x,  -  x  )2_S  ,  (3) 

&2  *  0 

where  u  and  u  are  the  surface  speeds  from  the  actual  distribution. 

el  e2 
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The  second  assumption,  based  on  physical  grounds,  is  that  the  momentum 
thickness  given  by  the  piecewise  wedge  flow  solutions  must  be  continuous.  The 
wedge  flow  relation  for  the  actual  momentum  thickness  9  is 


9  Re 


1/2  _ 


IX-X 

1  0 

(2  -  g) 

1/2 


(4) 


where  Re  is  the  Reynolds  number  based  on  body  axial  length,  0  is  the 
self -similar  momentum  thickness  defined  by. 


0  =  /  (1  -  f’)f'cn  .  (5) 

o 


f  is  the  solution  of  the  Falkner-Skan  differential  equation  and  q  is  the 
wedge  flow  similarity  variable  defined  by 


h 


Reu 

e 

(x-x  )(2-e) 

O 


1/2 


y 


(6) 


The  prime  in  Eq .  (5)  denotes  differentiation  with  respect  to  q. 

Equations  (2)  and  (3)  are  two  relations  for  the  determination  of  the 
three  constants  xQ,  k  and  g.  The  third  relation  is  obtained  by  evaluating 
Eq.  (4)  at  Point  1,  viz, 


Re 


1/2 


r  1 

x  ^  -x  ! 


—  !(2-S) 


“l  > 


1/2 


(7) 


and  is  known  from  calculations  on  the  preceding  interval.  The  only 
quantity  needed  of  the  three  is  g.  By  eliminating  x0  and  k,  one  obtains 
the  following  transcendental  equation  for  g  (see  Ref.  3  for  the  details): 
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6 

[l  +  (2  -  g)02z]2-6  =  1  +  $ 


(8) 


where 


u  -u 
e2  ®1 


(9) 


Z  = 


X2~X1 

9l‘Re0, 


(10) 


Re.  =  Re  •  u  •  3. 
°l  el  1 


(ID 


The  parameters  <p  and  Z  are  known  from  the  body  surface  speed  distribution  and 
the  boundary- layer  calculations  from  the  preceding  step.  In  Eq .  (8),  0  is  a 
function  of  0  which  requires  the  Falkner-Skan  solutions  to  be  known  and  an 
interpolation  scheme  to  be  used  to  determine  3.  The  interpolation  procedure 
will  be  discussed  later. 

Two  methods  are  available  for  the  solution  of  F,a .  (8)  for  8.  The  first 
involves  an  iteration  formula  obtained  by  taking  the  natural  logarithm  of 
both  sides  of  Eq .  (8)  and  solving  for  8  outside  the  logarithm.  This  formula 
is: 


2 

a+b  ’ 


a 


(12) 


where 
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a  =  Jln(l  +  <ji)  ,  (13) 

b  =  S.n[l  +  (2  -  B)02Z]  .  (14) 

Equation  (12)  has  excellent  convergence  properties  for  6  >  0,  converging  to 
four  decimal  places  in  S  in  four  to  six  iterations.  For  B  <  0,  successive 
approximations  to  S  oscillate  but  still  converge,  only  at  a  slower  rate  than 
for  B  >  0. 

The  second  method  is  a  straightforward  Newton  procedure  where  the  root 
of  the  following  equation  is  sought: 


f(b)  =  ~  *n[i  +  (2  -  e)e2z]  -  and  +  <j>)  =  0 


(15) 


Then  successive  iterates  for  B  are  given  by 


g(n+l)  =  g(n )  _ 


F(B) 


f'(b) 


(n) 


(16) 


In  Eq .  (16)  superscript  n  denotes  the  iteration  number  and  F  is  the 
derivative  of  F  with  respect  to  0  which  is  most  easily  determined  numerically. 
Equation  (15)  has  been  found  to  converge  in  all  cases  tried  and  in  fewer 
iterations  than  Eq .  (12). 

With  B  known,  @2  must  be  computed  to  continue  the  boundary-layer 
calculation.  The  formula  is  obtained  from  Eq.  (4)  with  the  help  of  Eqs.  (2), 
(3)  and  (7)  to  eliminate  k  and  xQ. 


The  result  is 
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9 


2 


1/2 

[i  +  (2  -  e)e2z] L  et 


(17) 


Note  that  over  the  Interval  xj  <  x  K,  x?  ,  the  Falkner-Skan  quantities  9 

A  a  a  A 

and  9  (and  hence  6*  and  H)  are  constant .  Thus  8,  9,  6  and  H  have  a  stair- 
step  behavior.  Consequently  the  physical  5  and  H  also  have  a  discontinuous 
behavior  with  only  9  being  smooth  because  it  was  forced  to  be. 


2.2  Extension  to  Suction 

The  extension  of  Smith's  method  to  include  suction  is  straightforward  and 
begins  with  the  wedge  flow  formula  for  vw,  the  velocity  normal  to  the  surface 
(c.f .  Ref.  1 ,  p.  235)  , 


Re 


(x-x  )(2-B) 
o 


1/2 


(18) 


where  fw  is  the  similarity  stream  function  evaluated  at  the  wall. 
Equation  (4)  can  be  used  to  eliminate  xQ  in  Eq.  (18)  giving  the  more 
useful  form: 


9  f 

w 

v  *  —  - 

w  9  Re 


(19) 


Because  9  and  f  are  constants  in  an  interval,  v 
w  w 

interval  junctions  but  this  behavior  will  not  turn 


will  be  discontinuous  at 
out  to  be  a  hindrance. 


1 
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In  the  present  application  just  the  right  amount  of  suction  must  be 
applied  so  that  small  disturbances  in  the  laminar  boundary  layer  will  not 
grow.  The  neutral  stability  curve  then  becomes  the  constraint  that  determines 
vw.  As  shown  by  Wazzan  et  al  (4] ,  neutral  stability  calculations  for  a  wide 
range  of  boundary-layer  flows,  when  plotted  as  (Reg*)crit  versus  the  shape 
factor  H,  form  a  single  curve  —  see  their  Fig.  1.  This  curve  in  tabular  form 
is  used  as  the  constraint  relation. 

As  illustrated  in  Fig.  1,  for  a  given  body  the  path  of  boundary-layer 
development  (Reg*  versus  H)  with  no  suction  proceeds  upward  toward  the 
critical  curve  as  the  distance  from  the  stagnation  point  increases.  The  point 
of  intersection  of  the  path  with  the  critical  curve  is  where  suction  must 
begin.  From  the  intersection  onward  the  path  is  constrained  to  follow  the 
critical  curve  and  thus  determines  the  minimum  suction  distribution  to 
maintain  laminar  flow.  Once  suction  is  turned  off  the  path  immediately  leaves 
the  critical  curve  and  transition  or  laminar  separation  or  both  follows 
quickly . 

On  each  interval  on  the  body  beginning  at  the  nose  the  Reynolds  number 
based  on  6*  is  calculated  and  tested  against  ( Re 5* ) c r i t  determined  from 
Wazzan's  relation  for  the  given  H.  Then  if 

(Rer*)  _  >  (Rer*)  .  _  , 

0  calc  0  crit 

the  constraint  relation  is  invoked  to  determine  6  and  fw.  This  relation  can 
be  written  as 


(Reg*),.  =  W(H) 


calc 


(20) 


where  W(H)  is  the  Wazzan  et  al  critical  Reynolds  number  relation.  Since  the 
calculated  boundary -layer  quantities  have  been  found  to  show  insignificant 
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differences  with  regard  to  where  Rej*  is  computed  in  each  interval,  provided 
the  interval  is  small,  the  beginning  point  of  the  interval  is  chosen  because 
the  calculation  of  8  is  simplified. 

Equation  (20)  is  therefore  written  as 

(Re6*)1=W(H)  .  (21) 

The  displacement  thickness  Reynolds  number  is  given  by 

( Re  <t* ) ,  =  Re  •  u  •  6*!  =  Re  •  u  •  9,  •  H 
6  1  ei  1  el  1 

(Re  =  Reg  •  H 

Then  Eq .  (21)  can  be  rewritten  as 

Re.  =  G(H)  ,  (22) 

91 

where 

G(H)  =  .  (23) 

n 

Equation  (22)  shows  why  the  beginning  of  the  interval  was  chosen.  Since 

Re  is  known  (because  9  is  continuous  and  hence  so  is  Re  ) ,  Eq .  (22)  provides 
91  9 
a  means  for  the  direct  determination  of  H  by  interpolation  of  the  relation 

of  Wazzan  et  al.  Thus  the  iteration  for  8  is  performed  along  an  H  =  constant 

line  which  eliminates  fw  as  an  independent  parameter,  as  will  be  seen.  The 

equation  that  determines  8  is,  of  course,  still  Eq.  (8). 
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Once  8  Is  known  for  the  interval,  62  is  computed  from  Eq .  (17)  and  then 

v  and  v  are  computed  from  Eq .  (19).  As  illustrated  in  Fig.  2,  the  typical 
W1  W2 

vw  distribution,  as  given  by  Eq.  (19),  has  a  saw-tooth  character.  A 
continuous  suction  distribution  is  usually  desired  for  design  purposes  as  well 
as  for  finite  difference  calculations  to  verify  the  approximate  results.  The 
best  agreement  between  approximate  and  finite  difference  boundary-layer 
calculations  was  found  for  a  "continuous"  distribution  of  vw  versus  x 
consisting  of  (vw)m^n  at  the  segment  junction  points,  i.e.,  the  largest 
suction  speed  magnitude. 

2.3  Axially  Symmetric  Formulation 

The  preceding  derivation,  as  Smith  shows,  can  readily  be  extended  to  the 
axially  symmetric  case  by  use  of  the  Mangier  transformation.  The  form  used  by 
Smith  is  as  follows: 


u  =  u 
e  e 


x  ir 


.2. 

°  I1 ! 


x  =  f  I  — •  dx 


y  =  r  y 


6=6 


(24) 


— . .  •  --- 


/ 


■ . . . 
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where  barred  quantities  denote  the  transformed  equivalent  two-dimensional 
flow,  y  is  distance  normal  to  the  body,  rQ  is  body  radius  and  R  is  maximum 
body  radius.  Then  in  the  wedge  flow  relations  all  physical  variables  are 
replaced  by  their  barred  equivalents  which  follows  from  the  fact  that  the 
axially  symmetric  boundary-layer  equations  transform  into  a  two-dimensional 
form  via  the  Mangier  transformation. 

By  definition,  the  momentum  area  is, 


6  =  / 
o 


1  - 


u 


r 

iT  dv 
u  R 

e 


(25) 


and  the  2-D  momentum  thickness  is, 


e  =  /  i 

o 


u 

u 

o 


J 


(26) 


The  two  are  therefore  related  as  follows: 


(27) 


Similarly,  the  displacement  area  is  related  to  the  displacement  thickness  by, 


5*  •  r  s*  • 


(28) 


Then  in  terms  of  9  the  critical  Reynolds  number  relation  becomes 


Re0, 


=  G(H)  , 


(29) 


/ 
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and  the  expression  for  suction  velocity  reads 


v 

w 


e 

e  Re 


(30) 


2 .4  Starting  Solution 

In  the  vicinity  of  the  stagnation  point  on  a  blunt  body,  the  relation  for 

9  is 


62  Re  =  3-  (2  -  B) 6“  .  (31) 

u 

e 


the  stagnation  point  x  =  0  and  =  0;  hence  Eq  •  (31)  is  indeterminate. 
Thus  6  at  the  stagnation  is  obtained  as  the  limit  of  Eq .  (31)  as  x  +  0 ,  viz. 


0.4283  I  Re 


du 

a*  j. 


^-1/2 


(32) 


where  the  values  used  for  8  and  9  at  the  stagnation  point  are: 


2  =  0.5  ,  0  =  0.3499  . 

so  sp 

Equation  (31)  must  also  be  used  to  determine  9  at  the  first  station 
downstream  of  the  stagnation  point,  x  =  x0 ,  because  the  general  procedure 
using  Eqs.  (2),  (3)  and  (7)  is  invalid. 


Thus 
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0 


2 


0.4285 


(33) 


2 . 5  Au tomat ion  of  the  Method 

The  original  Smith  method  was  set  up  for  graphical  solution  using  plots 
of  the  Falkner-Skan  solution  (fw  =  0)  in  the  form  of  A  versus  8  for  constant 
values  of  Z,  where  A  =  4>/Z.  For  solution  bv  digital  computer,  the  parameter 
$  is  more  convenient  and  the  graphs  must  be  replaced  by  a  large  data  base  of 
Falkner-Skan  solutions.  The  parameters  chosen  for  the  data  base  were  H  and 
6*  versus  8  and  fw,  including  the  separation  values.  The  ranges  covered  in 
8  and  fw  are  -  2.0  <  8  <  2.0  and  0  <  fw  <  2.5  which  has  been  found  to  be 
adequate  for  almost  all  the  numerical  examples  considered  to  date. 

The  curves  of  H  and  6*  versus  6  for  constant  values  of  fw  are  shown  in 
Figs.  3  and  4.  Each  of  these  plots  is  represented  in  the  data  base  by  227 
discrete  values.  The  data  base  also  contains  the  Wazzan  et  al  critical 
Reynolds  number  curve  represented  by  11  points  for  8  in  the  range  2.0  to 
3.0. 

To  obtain  the  values  of  H  and  6*  for  the  data  base,  the  Falkner-Skan 
differential  equation  was  solved  numerically,  for  8  and  fw  given,  using  the 
method  of  Keller  [5].  Double  precision  arithmetric  was  used  and  each  profile 
was  represented  by  50  intervals  in  the  similarity  variable  n  spaced  in  a 
geometric  progression. 

The  upper  limit  8  =  2  in  the  data  base  corresponds  to  the  maximum  value 
possible  for  wedge  flow  solutions  for  which  the  wedge  half  angle  is  180 
degrees.  Above  8=2,  the  present  method  breaks  down.  In  certain 


/ 
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applications,  notablv  for  flat  nose  bodies,  B  will  exceed  2  but  the  region 
in  x  over  which  this  occurs  is  usually  small.  Thus  setting  3  =  2  in 
such  a  region  produces  a  negligible  error  in  the  present  method,  as  noted 
by  Smith. 

The  condition  p  >  2  will  always  occur  close  to  the  nose  before  suction  is 
applied.  The  occurrence  of  3  =  2  can  be  detected  by  a  relation  derived  from 
Eq .  (8).  If  Eq .  (3)  is  written  in  the  following  form: 


*n[l+9"Z»(2-B)]  1  .  ,,  .  x 

- — -  =  -  in  (1  +  *>) 


(34) 


then  the  left-hand  side  is  seen  to  be  indeterminate  when  3=2.  By  letting 
e  -  2  -  B,  then 


£im 

E-0 


£n( l+4"Ze) 

€ 


"2 

e  z 


and  noting  that  for  3=2  and  fw  =  0,  9  =  0.2279,  the  final  limiting  form  of 
Eq .  (34)  is  found  to  be, 


Z  =  9 .6246  £n( 1  +  $) 


For  B  <  2,  the  following  inequality  must  therefore  hold: 

Z  >  9.6246  £n(l  +  4>)  .  (35) 

Thus,  whenever  this  inequality  is  violated  B  is  set  to  equal  to  2. 
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For  a  particular  approximation  to  R  in  the  solution  of  Eq .  (8),  for  given 
Z  and  p,  0  is  obtained  from  interpolated  values  of  H  and  6*.  These  latter 
two  quantities  are  determined  by  quadratic  Lagrange  interpolation.  For  the 
no  suction  case  (fw  =  0)  the  tables  of  H  and  5*  are  merely  interpolated  for 
the  given  R. 

For  the  case  fw  >  0,  H  is  first  determined  bv  linear  interpolation  of  the 
Wazzan  et  al  critical  Reynolds  number  curve  using  logjg  G.  The  values  of 
log^g  G  versus  H  are  given  in  Table  1. 


H 

log1Q  G 

2.0 

4.362 

2.1 

4.093 

2.2 

3.814 

2.3 

3.478 

2.4 

3.110 

2.5 

2.702 

2.8 

2.305 

2.7 

2.037 

2.8 

1 .863 

2.9 

1  .716 

3.0 

1.623 

Table  1 .  Critical  Reynolds  Number  Values  of  Vazzan  et  al. 


The  interpolation  procedure  to  obtain  6*  and  fw  is  as  follows: 

(1)  Interpolate  each  fw  =  constant  curve  in  the  data  base  to  obtain 
values  of  f  and  6*  for  the  given  H.  This  process  produces  curves  of 

A  JL 

6  versus  6  and  fw  versus  fi  for  H  =  constant. 

(2)  For  the  given  B,  interpolate  the  curves  from  Step  1  to  obtain 
6*  and  fw. 
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We  note  that  fM  is  actually  not  needed  until  the  6  calculation  has 
converged . 

The  Falkner-Skan  curves  of  H  versus  B  and  6*  versus  6,  with 
fw  =  constant,  in  the  region  8  <  0  and  fw  >  0  end  at  the  separation  point 
with  an  infinite  slope.  In  this  region  interpolation  along  H  =  constant  lines 
is  actually  a  much  more  accurate  way  to  determine  fw  and  6*  (for  a  given  8) 
than  interpolating  on  8  =  constant  lines. 

The  determination  of  8  and  fw,  as  dictated  by  the  critical  Reynolds 
number  curve,  will  automatically  preclude  the  occurrence  of  separation.  The 
reason  is  that  as  the  critical  Reynolds  number  increases,  H  must  decrease  and 

o 

will  always  be  less  than  the  corresponding  separation  value. 

If  suction  is  turned  iff,  H  will  increase  very  rapidly  and  transition 
will  ensue.  To  detect  this  occurrence  the  correlation  of  Wazzan  et  al  for 
Rex  versus  H  is  used.  The  Reynolds  number  is  based  on  the  body  arc  length,  as 
measured  from  the  stagnation  point,  of  the  point  where  a  disturbance 
amplification  ratio  of  e^  occurs.  This  correlation  has  been  found  by  Wazzan 
et  al  to  be  fitted  closely  by  the  equation: 

loglf)  [Rex(e9)]  =  -  40.4557  +  64.8066H  -  26.7538H2  +  3.3819H3  .  (36) 

If  Rex  at  the  particular  point  on  the  body  is  found  to  be  greater  than 
Rex(e^)  from  Eq .  (36),  then  transition  has  taken  place  and  the  boundary-layer 
calculation  is  stopped. 


i 
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3 .  MIMERICAL  RESULTS 


Suction  distributions,  as  well  as  boundary  layer  integral  properties, 
have  been  computed  for  four  different  bodv  shapes  at  Reynolds  numbers  ranging 

C  O 

from  10  to  10°.  Table  2  gives  the  particular  Reynolds  numbers  for  each  case. 


Case 

Bodv  Shape 

Reynolds  Numbers 

l 

Laminar  Flow  Bodv  in  Tunnel 

o 

o 

2 

NSRDC  Body 

106,  in7 

3 

Elat  Nose  Bodv 

106,  107,  108 

4 

F-57  3odv 

nC 

o 

X 

CM 

Table  2.  Bodv  Shapes  and  Reynolds  Numbers. 

These  particular  bodies  were  chosen  because  each  pressure  distribution  is 
significantly  different.  Case  1  has  no  adverse  pressure  gradient  near  the 
nose.  Case  2  has  a  nild  adverse  pressure  gradient  near  the  nose,  Case  3 
has  a  verv  strong  adverse  pressure  gradient  near  the  nose  and  Case  4  has 
a  medium  strength  adverse  pressure  gradient  near  its  mid-section.  In  both 
Cases  3  and  4  the  adverse  pressure  gradient  is  strong  enough  to  cause 
separation  of  the  laminar  boundary  layer  without  suction.  These  two  cases 
provide  the  most  severe  test  of  the  approximate  theory. 

To  verify  the  approximate  theory,  the  "smoothed"  suction  distribution, 
described  in  Section  2.2,  was  input  to  a  finite  difference  code  based  on 
Keller's  Box  Method  and  the  laminar  boundary-layer  properties  calculated. 
The  finite  difference  results,  which  can  be  considered  exact,  were  then 
compared  with  their  approximate  counterparts. 


i 
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Numerical  results  for  the  approximate  theory  were  obtained  on  a 
VAX  11/782  computer  using  single  precision  arithmetric.  Typical  CPU  times 
ranged  between  3  and  4  seconds  whereas  the  corresponding  finite  difference 
solutions  took  between  27  and  35  seconds  in  CPU,  a  factor  of  9  increase. 

The  step  size  distributions  along  the  bodies  studied  were  governed  by  the 
accuracy  requirements  of  both  the  inviscid  solutions,  from  which  the  pressure 
distributions  were  obtained,  and  the  finite  difference  boundary-layer 
solutions.  The  number  of  points  along  the  bodies  used  in  the  boundary-layer 
calculations  (approximate  and  exact)  ranged  from  73  to  88  over  a  total  arc 
length  range  of  0.63  to  0.72. 

Body  geometry  and  inviscid  surface  pressure  distribution  are  given  for 
each  case.  Then  comparisons  of  approximate  results  are  given  for 
distributions  of  shape  factor  H  and  momentum  area  9.  For  the  finite 
difference  results,  the  path  of  development  of  the  boundary  layer  is  given  in 
terms  of  log^g  (Re^*)  versus  H,  and  finally,  the  suction  velocity  distribution 
from  the  approximate  theory  is  given.  The  most  significant  of  the  above  as 
a  means  of  judging  the  worth  (accuracy)  of  the  approximate  theory  is  the  path 
of  boundary-layer  development.  The  closeness  with  which  the  finite  difference 
results  follow  the  critical  Reynolds  number  curve  as  a  result  of  suction  is  a 
measure  of  the  accuracy  of  the  approximate  theory  to  predict  the  correct 
suction  distribution. 

Results  for  the  laminar  flow  body  in  tunnel  are  shown  in  Figs.  5  through 
14.  The  body  geometry  appears  in  Fig.  5  and  the  pressure  distribution  in 
Fig.  6.  As  shown  in  Figs.  7  and  8,  the  agreement  between  approximate  and 
finite  difference  results  for  the  shape  factor  H  is  excellent  while  agreement 
for  momentum  area  9,  seen  in  Figs.  9  and  10,  is  good  with  a  maximum  difference 
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or  2.2  peruM  t  occurring  at  x  =  0.63  for  Re  =  10^.  The  finite  difference 
results  for  the  boundarv-laver  r>ath,  log  (Rej*)  versus  H,  closely  follow 
the  critical  Revnolds  number  curve,  although  slightly  above  it.  As  is  to 
he  expected,  increasing  the  Revnolds  number  moves  the  point  of  initial 
suction  application  forward,  from  x  =  0.495  at  Re  =  10^  to  x  =  0.045  at 
Re  =  10'.  This  dramatic  shift  is  a  consequence  of  the  streamlined  nose 
shape  and  favorable  pressure  gradient  induced  by  the  tunnel  blockage. 

The  N'SRDC  bodv  geometry  is  shown  in  Fig.  15,  its  pressure  distribution  in 
Fig.  16  and  numerical  results  in  Figs.  17  through  24.  The  streamlined  nose, 
which  fairs  into  a  cylindrical  midsection,  causes  a  mild  adverse  pressure 
gradient  to  exist  in  the  iuncture  region.  Again  the  agreement  in  H  at  both 
Revnolds  numbers  is  good  except  for  the  Re  =  10^  case  for  x  >  0.58  where 
suction  is  not  needed.  In  this  region  the  approximate  theory  predicts  a 
sudden  dip  in  H  which  is  less  pronounced  in  the  finite  difference  result. 

The  maximum  difference  in  H  between  the  two  in  this  region  is  4.3  percent. 

When  the  Re  =  10^  solution  was  first  run,  the  critical  Reynolds  number 
constraint  was  followed  for  x  >  0.58  which  resulted  in  blowing  rather  than 
suction  in  this  region.  A  test  was  then  put  into  the  code  to  prevent 
blowing.  Once  the  blowing  was  turned  off  for  x  >  0.58,  the  H  distribution 
was  found  to  undergo  the  sudden  drop  shown  in  Fig.  17.  The  comparison  in 
Fig.  17  thus  indicates  that  when  the  boundary  layer  undergoes  a  sudden  change 
the  approximate  theory  performs  poorly  at  first  but  soon  recovers.  A  further 
observation  about  the  behavior  of  the  H  distribution,  as  shown  in  Figs.  17 
and  18,  is  that  when  suction  begins  H  slowly  diminishes  as  x  increases.  This 
behavior  is  a  reflection  of  the  character  of  the  critical  Reynolds  number 
curve.  The  8  comparisons  for  this  case,  Figs.  19  and  20,  are  uniformly 
excellent  and  the  finite  difference  results  for  the  boundary-layer  path,  in 
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Figs.  21  and  22,  are  seen  to  follow  the  critical  Reynolds  number  curve 
reasonably  well,  although  slightly  above  it  as  in  Case  1.  The  "choppiness" 
of  the  finite  difference  results  is  attributed  to  a  numerically  noisy  pressure 
distribution.  Finally,  the  peaks  in  the  suction  velocity  distributions,  in 
Figs.  23  and  24,  correspond  to  the  location  of  the  maximum  adverse  pressure 
gradient  which  occurs  at  an  x  of  about  0.175. 

For  the  flat  nose  body,  the  geometry  appears  in  Fig.  25  and  the  pressure 
distribution  in  Fig.  26.  The  numerical  results  for  this  case  are  shown  in 
Figs.  27  through  38.  The  most  significant  feature  of  the  pressure 
distribution  is  seen  to  be  the  very  strong  adverse  pressure  gradient  present 
in  the  juncture  region  between  the  curved  part  of  the  nose  and  the  cylindrical 
raid  section.  Figures  27  through  29  show  the  agreement  for  H  to  be  good  at  all 
Reynolds  numbers  except  in  two  regions  on  the  body .  The  first  is  at  x  ~  0.03 
corresponding  to  where  3  >  2  in  the  finite  difference  solution.  As  already 
explained,  3  in  the  approximate  solution  cannot  exceed  2,  so  that  in  such  a 
region  it  is  held  constant  at  2.  Fortunately,  the  extent  in  arc  length  over 
which  this  occurs  is  very  small  (Ax  ~  0.02)  so  that  the  effect  on  the  solution 
downstream  is  negligible.  The  second  region  of  disagreement  is  in  the 
vicinity  of  the  maximum  adverse  pressure  gradient,  at  x  ~  0.1.  In  the  suction 
velocity  distributions,  shown  in  Figs.  36  through  38,  the  approximate 
theory  predicts  too  high  a  suction  peak  in  response  to  the  maximum  adverse 
pressure  gradient.  This  high  suction  peak  then  forces  H  in  the  finite 
difference  solution  to  dip  below  the  approximate  solution  -  see  Figs.  27 
through  29.  The  maximum  difference  in  H  between  the  two  is  5.7%  for 
Re  =  10^.  In  spite  of  the  departure  of  H  in  the  two  solutions,  the  agreement 
for  9  is  excellent  at  all  Reynolds  numbers,  as  shown  in  Figs.  30  through  32. 
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Finally,  the  boundary-layer  path.  Figs.  33  through  35,  clearlv  indicates  that 
the  suction  peak  at  x  ~  0.1,  from  the  approximate  theory,  is  too  high,  because 
the  finite  difference  results  suffer  an  excursion  on  the  stable  side  of  the 
critical  Revnolds  number  curve.  Afterward  they  follow  the  critical  Reynolds 
number  curve  reasonably  well,  although  slightly  on  the  high  side. 

The  final  case  is  the  F-57  low  drag  body  extensive lv  studied  by  Patel 
and  Lee  [6].  The  Revnolds  number  of  1 .2  x  10^  corresponds  to  the  value  in 
their  experiments.  The  body  geometry  is  shown  in  Fig.  39  and  the  pressure 
distribution  in  Fig.  40.  Patel  and  Lee  point  out  that  transition  on  this  bodv 
occurs  as  a  result  of  laminar  separation  followed  bv  turbulent  reattachment, 
in  the  vicinitv  of  x  ~  0.495.  .As  seen  in  Fig.  41,  the  approximate  theory 
predicts  a  raoid  rise  in  H  downstream  of  the  maximum  bodv  diameter,  in 
response  to  the  adverse  pressure  gradient.  Once  suction  begins,  at  x  =  0.46 
(slightlv  upstream  of  maximum  diameter)  H  begins  to  decrease  with  increasing  x 
as  a  result  of  following  the  critical  Reynols  number  curve.  On  the  other 
hand,  H  from  the  finite  difference  solution  of  x  >  0.45  rises  at  a  slower 
rate  but  finallv  coincides  again  with  the  approximate  solution  at  x  ~  0.6. 

As  in  the  flat  nose  case,  the  approximate  theory  does  not  perform  well  in  a 
region  where  the  pressure  gradient  changes  rapidly  from  favorable  to  adverse. 
The  9  distributions,  nevertheless,  agree  extremely  well,  as  Fig.  42  shows. 

The  boundarv-laver  path,  seen  in  Fig.  43,  reflects  the  behavior  seen  in 
Fig.  41.  The  point  where  suction  begins  in  the  finite  difference  results 
is  well  below  the  critical  Reynolds  number  curve  (in  the  stable  region).  This 
means  that  the  approximate  theory  predicts  too  much  suction  initially.  The 
initial  rise  in  suction  speed  in  Fig.  44  is  thus  too  rapid.  The  reason  whv 
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this  solution  was  not  run  any  further  rearward  on  the  bodv  than  x  =  0.66 
was  that  fw  at  that  point  had  nearly  exceeded  the  limit  of  the  data  base 
(2.5). 


4.  CONCLUSIONS 


The  extension  of  Smith's  piecewise  similarity  method  to  the  prediction  of 
the  minimum  distributed  suction  required  to  maintain  laminar  flow  in  an 
axisymmetric  boundary  layer  provides  an  automated  design  tool  that  is  both 
fast  and  reasonably  accurate  for  a  variety  of  body  shapes  and  Reynolds 
numbers.  In  cases  which  possess  pressure  gradients  that  change  rapidly  from 
favorable  to  adverse,  especially  when  the  adverse  pressure  gradient  is  strong, 
the  method  overpredicts  the  suction  speed  in  the  vicinity  of  the  maximum 
adverse  pressure  gradient.  Thus,  it  errs  on  the  side  of  conservatism. 

The  most  sensitive  (and  most  meaningful)  measure  of  the  accuracy  of  the 
approximate  method  is  the  boundary-layer  path,  log  (Ret*)  versus  H,  as 
obtained  from  a  finite-difference  solution  which  uses  the  suction  velocity 
distribution  from  the  approximate  theory.  For  the  four  bodies  studied,  the 
finite  difference  solution  was  found  to  follow  the  critical  Reynolds  number 
curve  quite  well  indicating  that  the  correct  amount  of  suction  was  being 
predicted.  As  already  indicated,  any  excursions  of  the  finite-difference 
results  from  this  curve  were  into  the  stable  region.  Finally,  in  all  cases 
examined,  the  agreement  between  approximate  and  finite-difference  results 
for  the  momentum  area  distribution,  9  versus  x,  was  excellent. 
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Figure  1 


Schematic  of  Path  of  Boundary-Layer  Development, 
log  (Re-*)  vs.  11. 
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Figure  2.  Schematic  of  Saw-Tooth  Behavior  of  Suction  Velocity. 
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PRESSURE  GRADIENT  PARAMETER,  p 

Figure  4.  Wedge  Flow  Displacement  Thickness  from 
Falkner-Sknn  Kquat ion . 
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Figure  7.  Boundary-Layer  Shape  Factor  Distribution, 
Laminar  Flow  Body  in  Tunnel,  Re  =  10^. 
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Figure  10.  Boundary-Layer  Momentum  Area  Distribution,  Laminar  Flow 
Body  in  Tunnel,  Re  =  10^. 
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Figure  12 


SHAPE  FACTOR,  H 

Path  of  Boundary-Layer  Development,  log  (Re^*)  vs.  H, 
Laminar  Flow  Bodv  in  Tunnel,  Re  =  10'. 


Figure  13.  Suction  Velocity  Distribution,  Laminar  Flow  Body 
in  Tunnel,  Re  =  10^. 


Figure  17.  Boundary-Layer  Shape  Factor  Distribution, 
NSRDC  Body,  Re  =  106. 
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Figure  20.  Boundary-Layer  Momentum  Area  Distribution, 
NSUDC  Body,  Be  =  1()7. 
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Figure  21.  Path  of  Boundary-Layer  Development,  log  (Re^*)  vs.  H, 
NSRDC  Bodv ,  Re  =  10<>. 


Figure  26.  Pressure  Distribution ,  Flat  Nose  Body 
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Figure  28.  Boundary-Layer  Shape  Factor  Distribution, 
Flat  Hose  Body,  Ke  =  107. 
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Figure  29.  Boundary-Layer  Shape  Factor  Distribution 
Flat  Nose  Body,  Re  =  10  . 


Figure  30.  Boundary-Layer  Momentum  Area  Distribution 
Flat  Nose  Body,  Re  =  10^*. 
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Figure  32.  Boundary-Layer  Momentum  Area  Distribution 
Flat  Nose  Body,  He  =  10  . 
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Figure  33.  Path  of  Boundary-Layer  Development,  log  (Ref*)  vs.  H, 
Flat  Nose  Body,  Re  =  10^. 
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SHAPE  FACTOR,  H 

Figure  34.  Path  of  Boundary-Layer  Development,  log  (Re**)  vs.  H, 
Flat  Nose  Body,  Re  =  10~. 


Figure  37.  Suction  Velocity  Distribution,  Flat  Nose  Body,  Re 
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Figure  38.  Suction  Velocity  Distribution,  Flat  Nose  Body,  Re 


Figure  41.  Boundary-Layer  Shape  Factor  Distribution,  F-57  Body, 


Figure  42.  Boundary-Layer  Momentum  Area  Distribution,  K-57  Body, 
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Figure  43.  Path  of  Boundary-Layer  Development,  log  (Re  *)  vs.  H, 
F-57  Body,  Re  =  1.2  *  106. 
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